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Abstract 

The effect of the surface coupling J s on the dependency of the layering transition 
temperature Tl as a function of the thickness N , of a spin— 1/2 Ising film, is studied using 
the mean field theory. It is found that for J s greater than a critical value (J sc = 1.30), 
the layering transition temperature decreases when the film thickness N increases for 
any values of the surface magnetic field H s . While, for J s < J sc , the behaviour of the 
layering transition temperature Tl, as a function of N, depends strongly on the values 
of H s . Indeed, we show the existence of three distinct behaviours of Tl, as a function 
of the film thickness N, separated by two critical surface magnetic fields H sc \ and H SC 2, 
namely: (i) for H s < H sc \, Tl increases with N; (ii) for H s & < H s < H SC 2, Tl increases 
for small values of N, and decreases for large value ones; (iii) while for H s > H SC 2, Tl 
decreases with increasing the film thickness. Furthermore, depending on the values of J s , 
the wetting temperature T W (T W = Tl when N — > oo for a given material), can be greater 
or smaller than the layering transition temperature of a film of thickness N of the same 
material. 
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1 Introduction 



The wetting and layering transitions of magnetic Ising systems have been studied by several 
authors. A simple lattice gas model with layering transitions and critical points has been in- 
troduced and studied in the mean field approximation by de Oliveira and Griffiths [1], Pandit 
et al. [2], Nightingale et al. [3] and Ebner et al. [4-7]. A variation of phase diagrams with 
the strength of the substrate potential in lattice gas model for multi-layer adsorption has been 
studied, using Monte Carlo simulations, by Patrykiejew et al. [8], and Binder and Landau [9]. 
One type of transitions is the layering transitions, in which the thickness of a solid film increases 
discontinuously by one layer as the pressure is increased. Such transitions have been observed 
in a variety of systems including for example 4 He [10,11] and ethylene [12,13] adsorbed on 
graphite. 

Ebner [14] carried out Monte Carlo simulations of such a lattice gas model. Huse [15] applied 
renormalization group technique to this model. It allowed the study of the effects on an atomic 
scale or order disorder transitions in the adsorbed layers, which may have considerable influence 
on the layering transitions and tracing back macroscopic phenomena on inter-atomic potentials. 
Benyoussef and Ez-Zahraouy have studied the layering transitions of Ising model thin films us- 
ing a real space renormalization group [16], and transfer matrix methods [17]. As for the film 
system, which is finite in one direction, it has been established that its magnetic properties can 
differ greatly from those of the corresponding bulk [18-22]. Experimental results [23] showed 
that the critical temperature of a vanadium film depends on the film thickness and its critical 
behaviour is like that of the two-dimensional system rather than that of the three-dimensional 
bulk. Although the great number of theoretical works made in the field of wetting phenomena, 
especially the effect of the nature of the substrate potential on wetting and layering transitions 
[1,2,16,17,24,25], the relation between the surface coupling and the wetting temperature is not 
yet sufficiently investigated. However, Hong [26] has studied the effect of the surface couplings 
on the behaviour of the critical temperature film as a function of its thickness. 
Our aim in this paper is to study the effect of the surface coupling (different than the bulk 
one) on the behaviour of the layering transition and the wetting transition temperatures, as a 
function of the film thickness of a spin— 1/2 Ising film, using the mean field theory. 
However, depending on the value of the surface coupling J s , the layering transition tempera- 
ture T L can increase or decrease with the film thickness and reaches the wetting temperature 
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T W (T W = Tl when N — > oo) for a sufficiently thick film. Furthermore, it is found that there 
exists a critical value H sc of the surface magnetic field below which the wetting temperature is 
independent on the surface coupling J s , while for H s greater than H sc the wetting temperature 
depends on the J s values. Such results can be useful for experiments. This paper is organised 
as follows. Section 2 describes the model and the method. In section 3 we present results and 
discussions. 



2 Model and method 

We consider a film formed with N coupled ferromagnetic square layers k — 1, ...,N (k — 1 is 
considered as a surface) infinite in the x and y directions, in an external and surface magnetic 
fields, as it is illustrated by Fig. 1. 

It is well known that the spin— 1/2 variable values can be written as: Si = (h/2)<Ji where 
<7j = ±1. For simplicity we assume, hereafter, that Ti/2 = 1. Hence, the Hamiltonian governing 
the system is given by: 

H. — - X] Jh3 a i a j -^HiCJi (1) 

<i,j> i 

The first sum runs over all nearest-neighbour sites. In several previous works [16,17,24,25], the 
case of a constant coupling J it j = =constant has been studied. Here, we consider the model 
where the constant couplings are given by: 



Ji,j — \ 



J s for e surface 

(2) 

Jb elsewhere 



In the last sum of the Hamiltonian (1), the total magnetic field Hi is applied on each site "i", 
and distributed according to: 



Hi 



H + H s for % e surface 

(3) 

H + H s /k a elsewhere 



with a being a positive constant. In all the following, we will be limited to a constant surface 
magnetic field H s applied only on the surface (k — 1). This is the case for a — > oo. 
The results of the model for a finite a value exhibits similar topologies of the phase diagrams 
as it was outlined by Binder et al. [27,28] and in some of our earlier works [16,17,29]. 



Using the mean field theory, we introduce the effective Hamitonian H = J2i ki<Ji + J2i Hi<Ji, with 
ki = J2j Jij a ji the sum runs over the nearest-neighbour sites 'f of the site V. The magnetization 
of each site i is given by: = T T7exp{-pHo) ' wnere Tr means the trace performed over all spin 
configurations. Since the system is invariant by translation in both x and y directions, = m k 
for each site 'i' of the layer k. Hence, the magnetisation of a layer k can be written as: 



mi = tanh(/9(4J s mi + Jbm 2 + Hi)) for k=l 

m k = tanh(/3(4J b m fc + J b m k ^i + J b m k+1 + H k )) for k=2,...,N-l 
mjv = tanh(/3(4 J^m^ + JbrriN-i + Hn)) for k=N 

With the free boundary conditions: m = m N+ i = 0. 

The total free energy of the system can be written as follows 



N 

k=l 



where F k stands for: 




(4) 



(5) 



-^Log(2cosh(f3(X 1 + H ± ))) + \mi\i for k=l 
-±Log{2cosh{(3{\ k + H k ))) + \m k \ k for k=2,..,N 

where 

Ai = 4J s mi + J b rfi2 for k=l 

A fc = 4J b m k + J & m fc _i + J h m k+l for k=2,...,N 
and j3 = \/(k B T) with T being the absolute temperature and ks the Boltzmann constant. 



(6) 



(7) 



3 Results and discussion 

Reduced values of the parameters T, T L) T W) H, H s and J s are investigated in this work. 
However, for simplicity the syntax "reduced" will be cancelled despite the fact that we use the 
notations X/ J b (X = T,T L ,T W , H, H s , J s ) and give numerical values of these reduced parame- 
ters. 

The notation \ k O N ~ k will be used for a configuration where the first top k layers from the 
surface, are in the state 'up' and the remaining iV — k bottom layers are in the state 'down'. 
In particular, 1^ (rep. N ) will denote a system where all the layers are with positive mag- 
netisation (rep. negative magnetisation). There exists iV + 1 possible configurations for a film 
formed with N layers. However, the transition of the layer "k" is characterised by the change 
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of the sign of its magnetisation while the magnetisations of the remaining N — k layers of the 
film keep their initial sign i.e.: i k ~ 1 () N ~ k+1 <-> l k O N ~ k . Hence, for k — 1 the transition is called 
the surface transition. 

The ground state phase diagram (T = 0) do not depend on the surface coupling constant J s . 
Hence, in the following we will give results and phase diagrams for T ^ 0. We solve numerically 
the equations (4) and (6) in order to establish the phase diagrams of the system. Indeed, Fig. 2 
shows that the layering transition temperature T L /J b , defined as the temperature above which 
the layers (k ^ 1) of the film begin to transit layer-by-layer, which is usually greater than the 
"surface transition" temperature. Except the special case of a film formed with two layers, 
where the surface transition temperature is close to the layering transition temperature T L /J b . 
It is understood that the wetting temperature T w /J b coincides with the limit of T L /J b for a 
sufficiently thick film (N — > oo). 

A numerical study of the layering transition temperature as a function of the film thickens is 
illustrated by Figs. 3a and 36 for two surface magnetic field values and selected values of the 
surface coupling constant. It is shown that for a smaller value of the surface magnetic field, 
the wetting temperature keeps a constant value (T w /J b 5.1 for H s /J b = 0.1) at any surface 
coupling constant values, as it is shown in Fig. 3a. For a higher value of the surface magnetic 
field (see Fig. 36), the scenario is inverted and the layering transition temperature T L /J b (as 
well as the wetting temperature) depends strongly on the surface coupling values. Indeed, for 
a fixed value of the surface magnetic field H s / J b = 1.0, some wetting temperature values are: 
T w /J b m 5.5 for a small value of the surface coupling constant J s /J b = 0.01 and T w /J b 2.4 
for a higher value J s /J b = 1.8, as it is summarized in Fig. 36. 

Nevertheless, there exists three classes of the layering transition temperature behaviours, namely: 
The behaviour (i) is a situation where the layering transition temperature T L / J b increases with 
the film thickness and stabilises at a certain fixed value. This behaviour is always seen for small 
values of the surface coupling and any surface magnetic field value. The behaviour (ii), where 
Ti/J b increases until a certain film thickness above which it decreases exhibiting a maximum. 
This situation is found for medium values of both the surface coupling and the surface magnetic 
field. Finally the behaviour (in), where T L /J b decreases continuously and stabilises at a value 
close to T w / J b . This is generally the case for higher values of surface coupling. 
It is worth to note that the wetting temperature T w /J b (which coincides with T L /J b for a suffi- 
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ciently thick film) decreases as H s / Jb increases at a fixed surface coupling constant J s j Jf,. This 
is because the surface region becomes magnetically harder than the deeper layers, and as the 
film is thicker, the strength of the internal magnetic field decreases leading to a decrease of 
T w / Jb- 

We found two critical surface coupling constants J sc i/Jb and J SC 2/Jb- Jsci/Jb separates the 
behaviours (i) and (ii), while J SC 2/Jb separates the behaviours (ii) and (in). The dependency 
of these critical parameters, J sc i/Jb and J sc2 / 'Jb, as a function of surface magnetic field is plot- 
ted in Fig. 4. For small values of the surface magnetic field H s / J b , the constants J sc i/Jb and 
J sell 1 Jb become similar and are close to the value 1.30, and vanish as the surface magnetic field 
H s / J h increases. 

The results of the Figs. 3(a,b) are summarised in Fig. 4, since for small values of surface 
magnetic field H s /Jb the three behaviours of Tl/ Jb, namely (i), (ii) and (Hi) are found in this 
figure, when increasing the surface coupling J s /Jb values. 

The two critical surface coupling constants J sc i/Jb and J SC 2/Jb correspond, respectively, to two 
critical surface magnetic fields H scl /J b and H sc2 / Jb separating the three regimes of the T L /J b 
behaviours. For sufficiently large values of J s / Jb (J s /Jb > 1.30), Tl/ Jb decreases when increas- 
ing the film thickness N for any values of H s / J b . While, for J s j J b < 1.3, the layering transition 
temperature T L /J b exhibits the three behaviours: (i), (ii) and (iii) depending on the H s / J b 
value. In particular for H s /Jb > H sc2 /Jb, TL/Jb is always decreasing when the film thickness 
increases, see Figs. 3a, 3b. It is also interesting to examine the surface coupling effect on the 
layering transition temperature, for a fixed surface magnetic field value. As it is illustrated 
in Fig. 5, for a film size N = 5 layers and several surface magnetic field values, the layering 
transition temperature increases for small values of H s /Jb and decreases when H s /Jb reaches 
higher values. This figure shows also that the layering transition temperature decreases when 
increasing the surface magnetic field H s / J b , for a fixed value of the surface coupling J s j J b . This 
result was already outlined in earlier works e.g. Pandit et al. [2] and in some of our recent 
works [24,25]. 
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4 Conclusion 

In conclusion, we have studied the effect of the surface coupling J s on the wetting and 
layering transition temperatures as a function of the thickness N, of a spin— 1/2 Ising film, 
using mean field theory. We showed the existence of a critical value (J sc = 1.30) of the surface 
coupling constant J s , above which the layering transition temperature Tl decreases when the 
film thickness increases for any values of the surface magnetic field H s . While, for J s < J sc , 
there exists three distinct behaviours, namely: (i) for H s < H sc i, T L increases with N; (ii) 
for H sc \ < H s < H SC 2, Tl increases for thin films, and decreases for thick ones; (iii) while for 
H s > H sc2 , T L decreases with increasing the film thickness. Furthermore, for J s > J sc the 
wetting temperature T W (T W = T L when N — > oo), is independent on J s . While for J s < J sc T w 
depends on the value of J s . Such result may be useful for experiments. 
References 

[1] M. J. de Oliveira and R. B. Griffiths , Surf. Sci. 71, 687 (1978). 

[2] R. Pandit, M. Schick and M. Wortis, Phys. Rev. B 26, 8115 (1982); 
R. Pandit and M. Wortis, Phys. Rev. B 25, 3226 (1982). 

[3] M. P. Nightingale, W. F. Saam and M. Schick, Phys. Rev. B 30,3830 (1984). 

[4] C. Ebner, C. Rottman and M. Wortis, Phys. Rev. B 28,4186 (1983). 

[5] C. Ebner and W. F. Saam, Phys. Rev. Lett. 58,587 (1987). 

[6] C. Ebner and W. F. Saam, Phys. Rev. B 35,1822 (1987). 

[7] C. Ebner, W. F. Saam and A. K. Sen, Phys. Rev. B 32,1558 (1987). 

[8] A. Patrykiejew A., D. P. Landau and K. Binder, Surf. Sci. 238, 317 (1990). 

[9] K. Binder in Phase Transitions and Critical Phenomena, edited by C. Domb and J. L. 
Lebowits (Academic, New York, 1983), Vol. 8. 

[10] S. Ramesh and J. D. Maynard, Phys. Rev. Lett. 49,47 (1982). 

[11] S. Ramesh, Q. Zhang, G. Torso and J. D. Maynard, Phys. Rev. Lett. 52,2375 (1984). 



7 



[12] M. Sutton, S. G. J. Mochrie and R. J. Birgeneou, Phys. Rev. Lett. 51,407 (1983); 

S. G. J. Mochrie, M. Sutton, R. J. Birgeneou, D. E. Moncton and P. M. Horn, Phys. Rev. 
B 30,263 (1984). 

[13] S. K. Stija, L. Passel, J. Eckart, W. Ellenson and H. Patterson, Phys. Rev. Lett. 51,411 
(1983). 

[14] C. Ebner and W. F. Saam, Phys. Rev. A 22, 2776 (1980); ibid, Phys. Rev. A 23,1925 
(1981); ibid, Phys. Rev. B 28,2890 (1983). 

[15] D. A. Huse , Phys. Rev. B 30,1371 (1984). 

[16] A. Benyoussef and H. Ez-Zahraouy, Physica A, 206, 196 (1994). 

[17] A. Benyoussef and H. Ez-Zahraouy, J. Phys. I France 4, 393 (1994). 

[18] A. J. Freeman, J. Magn. Magn. Mater, 15-18, 1070 (1980). 

[19] R. Richter, J. G. Gay and J. R. Smith, Phys. Rev. Lett. 54, 2704 (1985). 

[20] R. H. Victora and L. M. Falicov, Phys. Rev. B 31, 7335 (1985). 

[21] C. Rau, C. Schneider, G. Xiang and K. Jamison, Phys. Rev. Lett. 57, 3221 (1986). 

[22] D. Pescia, G. Zampieri, G. L. Bona, R. F. Willis and F. Meier, Phys. Rev. Lett. 58, 9 
(1987). 

[23] C. Rau, G. Xiang and C. Liu, Phys. Rev. Lett. A 135, 227 (1989). 

[24] L. Bahmad, A. Benyoussef and H. Ez-Zahraouy, Phys. Rev. E 66, 056117 (2002). 

[25] L. Bahmad, A. Benyoussef and H. Ez-Zahraouy, M. J. Condensed Matter 4, 84 (2001). 

[26] Q. Hong, Phys. Rev. B 41, 9621 (1990); ibid Phys. Rev. B 46, 3207 (1992). 

[27] K. Binder and D.P Landau, Phys. Rev. B 46, 4844 (1992) 

[28] K. Binder in Cohesion and structure of Surfaces vol. 4, Elsevier (1995). 

[29] L. Bahmad, A. Benyoussef A. Boubekri and H. Ez-Zahraouy, Phys. Stat. Sol. b 215, 
1091 (1999). 

8 



Figure Captions 



Figure 1. A sketch of the geometry of the system formed with N layers and subject to 
the surface magnetic field H s and the global magnetic field H . J s denotes the interaction cou- 
pling constant between the spins of the surface; J b is the interaction coupling constant between 
the spins of the bulk as well as between the spins of the bulk and those of the surface. 

Figure 2.: Phase diagram, in the (H/ J b ,T/ J b ) plane, showing the layering transitions, and 
the definition of the layering transition temperature Ti/Jb, which is different from the surface 
transition temperature. This figure is plotted for a system size with N = 10 layers, a surface 
magnetic field H s /J b = 1.0 and a surface coupling constant J s /Jb = 1-5. 

Figure 3.: Layering transition temperature, T L / J b) behaviour as a function of the system size 
( number of layers) with a) H s / J b = 0.1 and b) H s / J b = 1.0. The number accompanying each 
curve denotes the surface coupling J s / J b value. 

Figure 4.: Critical surface couplings J sc i/J b and J sc2 / J b profiles as a function of the surface 
magnetic field H s /J b . The behaviours {%) and (ii) are separated by the J sc i/J b line; while 
J S c2/ Jb separates the behaviours (ii) and (in). 

Figure 5.: Layering transition temperature behaviour as a function of the surface coupling 
J s /J b for N = 5 layers and selected values of the surface magnetic field H s /J b : 0.1, = 0.2, 0.3, 
= 0.4, 0.5 and 0.6. 
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